function stress_max=post_read_stress(dt,n)
stress_max=[];
parfor i=1:n
    z=1;
    stress=[];
    f=fopen(['result\101_',num2str((i-1)*dt),'_',num2str(i*dt),'.pch']);
    while true
        tline1=fgetl(f);
        if length(strfind(tline1,'HEXA'))==1
            while true
                tline1=fgetl(f);
                if tline1(1)=='$'
                    break
                end
                for j=1:2
                    fgetl(f);
                end
                tline1=fgetl(f);
                stress(z)=sscanf(tline1(60:72),'%f');
                z=z+1;
                for k=1:8
                    for j=1:6
                        fgetl(f);
                    end
                    fgetl(f);
                end
                for j=1:4
                    fgetl(f);
                end
            end
            break
        end
    end
    fclose(f);
    stress_max(i)=max(stress);
    find(stress==stress_max(i));
end
delete(gcp('nocreate'));
end